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Abstract 

The global in-time semiclassical and relaxation limits of the bipolar quantum 
hydrodynamic model for semiconductors are investigated in . We prove that 
the unique strong solution converges globally in time to the strong solution 
of classical bipolar hydrodynamical equation in the process of semiclassical 
limit and to that of the classical Drift-Diffusion system under the combined 
relaxation and semiclassical limits. 
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1 Introduction 

Recently, the quantum hydrodynamic(QHD) model for semiconductors is derived and 
studied in the modelings and simulations of semiconductor devices (like MOSFET and 
RTD) in ultra-small size (say nano-size), where the effects of quantum mechanics, such as 
particle tunneling through potential barriers and built-up in quantum well, are taken into 
granted and dominate the transportation of electron and/or hole under the self-consistent 
electric field. 

The basic observation concerning the quantum hydrodynamics is that the energy den- 
sity consists of one additional new quantum correction term of the order 0{h) introduced 
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first by Wigner [32] in 1932, and that the stress tensor contains also an additional quantum 
correction part [2|, [3] related to the quantum Bohm potential (or internal self-potential) [Ij 

with observable p > the density, m mass, and h the Planck constant. The quantum 
potential Q was introduced by de Broglie and explored by Bohm to make a hidden vari- 
able theory and is responsible for producing the quantum behavior, so that all quantum 
features are related to its special properties. Such possible relation was also implied in the 
original idea initialized by Madelung [27] in 1927 to derive quantum fluid-type equations, 
in terms of Madelung's transformation applied to wave function of Schrodinger equation of 
pure state. In fact, based on this idea, one is able to derive quantum fluid type equations 
from the (nonlinear) Schrodinger equation of pure-state [101 [E] • 

The moment method is employed recently to derive quantum hydrodynamic equa- 
tions for semiconductor device at nano-size based on the Wigner-Boltzmann (or quantum 
Liouville) equation [29] 

Wt + ^- + ^F^W = [Wt]c (1.2) 
m 

where W = W{x,^,t), {x,C,t) e x x is the distribution function, and P the 
pseudo-differential operator defined by 

The electrostatic potential $ = ^{x,t) is self-consistent through Poisson equation 

XoA^ = q{ J Wd^-C), 

with Ao > the permitivity characteristic of device, q the elementary charge, and C = 
C{x) > the given doping profile [29], and [Wt]c refers to the collision term. In fact, 
applying moment method to the Wigner-Boltzmann equation (11.21) near the "momentum- 
shifted quantum Maxwellian" [32] together with appropriate closure assumption [H [11] , 
one can obtain the quantum hydrodynamic equation [B] . For more derivation and related 
topics on the modelling of quantum models, one refers to [2^ \TU\ [S] and the references 
therein. 

In the present paper, we consider the bipolar quantum hydrodynamic model of semi- 
conductors (for carriers of two type) 

dtPi + V ■ (piUi) = 0, (1.3) 

dtipiUi) + V ■ {piUi u^) + VP^{p^) = QiPiE + i-p^v(— ^) - (1.4) 

2 ^/Pi Ti 



AV ■ E = Pa- pb-C, V X E = 0, E{x) 0, |x| +00, (1.5) 
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where {x,t) E x i?+ and the index i = a,b and Qa = I, qt = —1. The observable 
Pa > 0,pb > 0,Ua,Ub and E are the densities, velocities and electric field, respectively. 
Pa{-),Pb{-) are the pressure- density functions. The parameters e > 0, Ta = u = r > 0, 
and A > are the scaled Planck constant, momentum relaxation time, and Debye length 
respectively. C = C{x) is doping profile. 

In the real simulations of semiconductor devices, the size of the device is rather small 
(in nanosize, for instance). This in turn makes the scaled parameters r, e, A rather smaller 
due to different situations under consideration [301 [29]. general, the scaled parameters 
£, r, A are expressed as 

2 _ ^ ■> 2 _ ^O^^J 2 _ ^bTqTq 

^ ~ 2mKBToL^ ' ~ Nq^L^ ' ^ ~ mL2 

where we recall that the physical parameters are the elementary charge q, the Boltzmann 
constant fc^, the elective electron mass m, the reduced Planck constant h, the permittiv- 
ity Ao, the ambient temperature Tq, and the characteristic device length L and density 
A^. The typical values of the parameters for semiconductors are given in [29]. There- 
fore, one of the both mathematically and physically important problems is to justify the 
asymptotic approximation (or behavior) of the macroscopic observable of the quantum 
hydrodynamical model subject to the small parameters mentioned above. 

In the present paper, we investigate the asymptotical analysis with respect to the 
scaled small parameters of bipolar time-dependent quantum hydrodynamical model. To 
begin with, let us present a complete description about the small-scale asymptotics of the 
QHD model. We first consider the semiclassical limit. Let e —>■ formally in fll.3p - fll.5l) . 
we get the well-known bipolar hydrodynamic (HD) model [HI d] 

dtPi + V ■ (piUi) = 0, (1.6) 

OIL' 

OMui) + V ■ (piUi ® Ui) + VPiipi) = qiPiE - (1.7) 
AV ■ E = Pa- pb-C, V X E = 0, E{x) 0, \x\ +00. (1.8) 

This limiting process shows the semiclassical approximation of bipolar quantum hydrody- 
namical model in terms of bipolar hydrodynamical model for small Planck constant, and 
describes the relation from quantum mechanics to the classical Newtonian mechanics. 

The semiclassical limits of the stationary unipolar quantum hydrodynamical model 
(carrier of one type) are well studied recently. In one dimensional bounded domain, the 
semiclassical limit of the thermal equilibrium solutions [5] and the isentropic subsonic so- 
lutions [12] are analyzed respectively due to different boundary conditions. This limit is 
also investigated for a stationary unipolar viscous quantum hydrodynamical system [7] for 
a special class of viscosity in one-dimensional interval subject to the boundary condition 
of density and quantum Fermi potential, where the communication between vanishing 
viscosity and semiclassical limit is also investigated in subsonic regime. For bipolar sta- 
tionary quantum hydrodynamical model, the semiclassical limits are investigated in multi- 
dimensional bounded domain for isothermal solutions in thermal equilibrium state |3n[23] . 
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by recovering the minimizer of limiting functional of a quantized energy functional corre- 
sponding the original system, and in multi-dimensional unbounded domain for stationary 
isentropic system [35]. A rigorous analysis is also made for the bipolar viscous quantum 
hydrodynamical system |23j . 

However, all those analysis for stationary problems can not apply to the time-dependent 
case because that unlike the case of stationary problem, the maximum principle usually 
does not apply to the time- dependent case and it is not clear how to derive enough the a- 
priori estimates with respect to time (derivatives) so as to pass into the semiclassical limit. 
Although such process of semiclassical limit has been investigated recently for nonlinear 
Schrodinger equation [2^1^ for potential flow in terms of Friedrich-Kato-Lax's theory and 
is concerned with the finite (short) time theory, the frame work does not apply here to 
general multi-dimensional rotational (non-potential) flow and is not fit for global in-time 
theory. We should do the semiclassical limit for QHD model in a different way in order 
to present the global in-time semiclassical limit for general non-potential flow. 

Next, we turn to the analysis of relaxation limit. To this end, let us introduce the 
diffusion scaling as [201 EH] 

t u- t 

^^^^ t^-^ {plulE-){x,t) = {p„^,E){x,-). (1.9) 

T T T 

Then (ll.Sp - lll.Sp can be rewritten as 

dtpl + V-ipluJ) = 0, (1.10) 

r'dtiplul) + ■ (p[< ® ul) + VP.(p[) = q^plE^ + f p[V(^) - p[<, (1.11) 

-E^ = pI- pI- C(x), V X = 0, E^{x) 0, |x| +oo, (1.12) 

Also formally, let r — > in (I1.10l) - fll.l2l) . the quantum drift-diffusion(QDD) model is 
obtained 

^tp^ + V{q,.p,E - VPiipi) + ^p,V(^)] = 0, (1.13) 

A^V ■ E = Pa- pb- C{x), V X E = 0, E{x) 0, \x\ +oo, (1.14) 

This limiting process provides a singular approximation of quantum hydrodynamical 
model via parabolic quantum Drift-Diffusion model for small momentum relaxation time. 
Note that although there are many results obtained for classical hydrodynamic model 
[H [281 [22], few is known for the relaxation limit for the quantum hydrodynamical 
model due to the less of enough information to control the nonlinear third order disper- 
sion term. Although the relaxation limit of the stationary solutions are investigated in 
one-dimensional bounded domain for unipolar case [12], and in multi-dimensional bounded 
domain for bipolar case [23], like the situation of semiclassical analysis, all these studies 
seems not enough in the resolution of the time-dependent problems. Note that, the singu- 
lar relaxation time limit presented above is not mathematically rigorous, the first rigorous 
analysis result about relaxation time limit of QHD model has been obtained recently in 
[20] . where the QHD system is proven to be approximated by a quantum Drift-Diffusion 
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model(QDD), a nonlinear parabolic equation, for small relaxation time. However, this 
analysis depends strongly on the effects of the nonlinear dispersion. That is, the scaled 
Planck constant is required to be fixed in order to help getting enough control to pass 
into the relaxation limit. This analysis is therefore not enough to prove the relaxation 
limit for possibly arbitrary small Planck constant e. Thus, it is natural for us to consider 
the relaxation limit of quantum hydrodynamical model for any small Planck constant e 
and furthermore the combined relaxation and semiclassical limit. In fact, we can show in 
the present paper that one can derive the following limiting drift-diffusion (DD) model 

5tp. + Vfep.E - VP.(p.)] = 0, (1.15) 

-E = Pa- pb- C{x), V X E = 0, E{x) -> 0, l^l +00, (1.16) 

by setting r ^ and e — in ( 11.10p -f fLT2l) for strong solutions. Note here that although 
we only deal with the combined relaxation and semiclassical limits for the quantum hy- 
drodynamical model f ll.l0p -( n..l2p . we claim that the analysis made here does not require 
any (communication) restriction between e and r. That is, one can fix any of the two 
parameters e and r and let the other tend to zero. 

We shall also mention the asymptotical analysis about the zero-Debye length limit for 
QHD model. This process is quite well understood for both stationary problems [31] [12] for 
one and multi-dimension bounded domain respectively and the time- dependent problem 
for mult i- dimension [25j. We omit the corresponding analysis here. 

The rest part of the paper is arranged as follows. The main results related to semiclas- 
sical limit and relaxation time limit are presented in section 2, the proofs are established 
in section 3. 

Notations: C or c always denote the generic positive constants. L^(M^) is the space of 
square integral functions on with the norm || ■ || or || ■ ||L2(]g3). H''{M.^) with integer 
k > 1 denotes the usual Sobolev space of function / satisfying d^f G L^(M^)(0 < i < k) 
with norm 

y 0<|a|<A: 

here and after a G N^,D'^ = d^^d^ld^!^^ for |a| = si -|- S2 + S3, Especially || • ||o=|| ■ ||. 
Let ;B be a Banach space, C'^([0, t]; B) denotes the space of ;B- valued k-times continuously 

differentiable functions on [0,t]. We can extend the above norm to the vector- valued 

3 

function u = {ui,U2, ,^3) with \D°'u\^ = ^ \D°'Ur\'^ and 

r=l 

II D'ur= [ (X^ 5^ (D^o'M^, 

r—i |Q!|=fe 

k 

and II u \\k=\\ u \\jjk(RS)= II D'u ||, || / ||L-'([o,r];i3)= sup || f(t) \\b. We also 

j=0 0<t<T 

use the space n''{M.^) = {/ G L%R^),Df G H''-\R^)},k > 1. Sometimes we use 
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II (., ., ...) ||i^fc(M3) or II (., ., ...) Ilfc to denote the norm of the space H^{R^) x H^{R^) x ■ ■ 
■ X H''{R^) and the n''{R^) as well. 

2 Main results and Preliminary 
2.1 Main results 

We consider the initial value problem for the quantum system fll.3l) - fll.5p with following 
initial data 

{pi,Ui){x,0) = {pio,Uio){x), pio{x) ^ p*, Ui,,{x) 0, \x\ ^ +00 (2.1) 

with i = a,b. From now on, we set the scaled Debye length to be one A = 1 for simplicity. 

First of all, we have the global existence and uniqueness theory of the IVP problem 
for the quantum system fll.3p - fll.5l) and (12. ip . 

Theorem 2.1 (Global existence) Let the parameters e > 0,t > be fixed. Assume Pa, 
Pb G C^(0,+cx3) and C{x) = c* is a constant satisfying for two positive constants pi, pi 
that 

pl-pl-c* = 0, PM), n'(P6)>0, (2.2) 
Suppose pao > 0,Pbo > and ./p^ - ^/f^, - ^/f^, Uao, Ubo) e i:/'^(M^) x7i^(M3). Then, 

there zs Ai > so that if Aq := \\{^/p^ - y/Pb^ - \^,Uao,Ubo)\\H^xnHm ^ ^i' 
unique solution (p^, pi, u^, ttf , E'^) of the IVP problem i\1.3h — (I i.jp and {\2.1\} exists globally 
in time with p^, p| > and satisfies 

{pl - pI, e') e c\o, T- H'^-^''(w^)), ul e c\o, T- n^-^\w^)), k = 0,1,2 

and 

II (Pa - Pl^Pb - pDh^m + ||^1|l-(IR3) + ||«,Mfe)|U«>(M3) — > 0, 

as time tends to infinity. 

Remark 2.2 Unlike the unipolar quantum hydrodynamical model fTR l2R [7^1 H^, we 

can not get the exponential convergence to the asymptotical equilibrium state for bipolar 
quantum model due to the coupling and cancelation of two carriers. Usually, the optimal 
decay rate is algebraic and is left for the further research I34j 

We then state semiclassical limit e 0+ of the global in time solutions to the 
IVP (ll.3p - (ll.5p and (12. ip for any fixed momentum relaxation time r > 0. 

Theorem 2.3 (Global semiclassical limit) Let r = 1 and {p^, pl,ul,ul,E^) be the 
solution of the IVP problem ( I i . gp -( TOl) and ^2.1\\ given by Theorem \2.1\ Then, there is 
{pa, Ua, Pb, Ub, E) with pa > 0, pb > SO that as the Planck constant e ^ 0, it holds 

pl^Pi strongly m C{0,T;Cl n H^-'); u\ ^ m strongly m C{0,T;Clv\n]-^'); 
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strongly m C(0, T; n H^c). s e (0, i). 

for any T > 0, i=a,h. Note here that {pi,Ui, E) with i = a,b is the global in-time solution 
of I VP problem of the bipolar hydrodynamic model (1 i.^p - fl j.<^ and {\2.1\i . 

Finally, we consider the combined semiclassical and relaxation limits for the quan- 
tum hydrodynamical model fll.3p - fll.5p . To this end, we consider indeed the initial value 
problem for the re-scaled system fll.10p - fll.12p together with the following initial data 

(p[,<)(x,0):=(p[„,O = (p„^)(x). (2.3) 

It is easy to verify that there is a unique global in-time strong solution {p^^'^\ u^i'^\ -E*^^'^-*) 
with i = a,b for the IVP problem f ll.l0p -f |1.12l) and (12.31) based on the Theorem 12.11 and 
the diffusion scaling (11. 9p . What left is to establish the uniform estimates with respect to 
the parameters e > 0, r > in order to pass into the limits. We have. 

Theorem 2.4 (Global relaxation and semiclassical limits) Let {pf'^\u'['^\E'^'^''^'>) 

with i = a,b be the unique global solution of the bipolar QHD equations f li.iOp - O.i^P and 
( \2. 51) obtained in Theorem \2.1\ then there exist {pa, Pb, E) such that as e — and r — > 

strongly tn C(0, T; n i/f-^(i?=')), 
^ E strongly m C7(0, T; Cl H Hf-'{R^)), 
rVf strongly m L\0, T; W^^^{R^)), s G (0, i). 

and {pa, Pb, E) is the strong solution of the IVP problem of bipolar Drift-Diffusion system 
(1 j. j5|) - fli.igp and initial data {pa, Pb){x,0) = (paOiPbo)- 

Remark 2.5 Although we only state the combined relaxation and semiclassical limits for 
the quantum hydrodynamical model i\1.10^ - i\1.12\) here, we claim that the analysis made 
here does not require any (communication) restriction between e and r. That is, one 
can fix any of the two parameters e and r and let the other tend to zero. Moreover, our 
analysis for the bipolar model (\1.10^ - (\1.1^ can be applied to justify the semiclassical limit 
and relaxation limit for the unipolar model 14 , 



Remark 2.6 Although we have only taken the steady state of constant solution in the 
profile in above theorems, we claim that our analysis below is valid for general subsonic 
steady state. 

2.2 Some lemmas 

Lemma 2.7 Let f G H^{R^), s > |. There is a unique solution of the divergence equation 
V ■ u = f, VxM = 0, u{x) 0, +oo. 

satisfying 

II U \\l^{R3)< C II / ||L2(i?3), II DU ||//.(i?3)< C II / ||H»(ij3) ■ 
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Lemma 2.8 Let f G H^{R'^), s > | with V ■ / = 0. There is a unique solution u of the 
vorticity equation 

V X u = f, V ■ u = 0, u{x) — * 0, |x| — > +00. 

satisfying 

II U ||l6(r3)< C II / ||l2(^3), II DU ||/fs(/j3)< C II / ||//s(ij3) . 

We will also use the Moser type calculus lemmas. 
Lemma 2.9 Let f,g e H'{R^)f]L'^{B:^), then it holds 

II D'^ifg) \\< C II g ■ \\ f \\ +C || / ■ \\ D^g \\ 

for a G iV^ 1 < |a| < s, s > is an integer. 

Lemma 2.10 Let f G H^{E?) with s > be an integer and function F[p) smooth enough 
and F(0) = then F{f){x) G H'{R^) and 

II Fif) \\h={r3)< C II / \\h^{r3) ■ 

3 The proof of main results 

The local in-time existence result of QHD model has been obtained in [T^ . The 
framework used there is to study an extended problem derived based on a deposition of 
the original problem, which in turn implies the expected problem as a special case. The 
method employed in [26l [T5] can be applied to our bipolar model directly. The proof is 
straightforward, and we have 

Lemma 3.1 (Local existence) Let the parameters e > 0, r > 0, A > be fixed. Assume 
that there are constants p*, > and c* satisfying pl — pl—c* = 0, C{x) — c* G H^{R^), 
and Pa, Pb e C^(0,+oo). Assume {^/pl^ - ^/f^,Uio) G H^{R^) x n^{R^) with > 0, 
then there exists a finite time T* > such that the unique solution {pa, Pb,Ua,Ui,, E) with 
Pa > 0, Pb > of the problem {\1.!!J\} - {\1.5\} and exists in [0,T*], and it satisfies 

Pi- p* eC^{[Q,T*]-H'^-^\R^)), UieC\[Q,T*]-n^~^^{R^)), A; = 0,1, 2. 

E eC^{\Q,T*]-7f-^{R^)), A; = 0,1. 

Here, we also mention the global existence theory for the quantum hydrodynamical 
model and bipolar hydrodynamical model. The well-posedness of steady state subsonic 
solutions has been proved also in [T9l [33] . Transient solutions are shown to exist either 
locally in time [I3l HI] or globally in time for data close to a steady state [I5l [161 [El 126] . 
The bipolar hydrodynamic(HD) model of the global solutions has been studied in [9]. 
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3.1 Reformulation of original problem 

In this section, we study the global solutions and the asymptotic limits with the case 
C{x) = c*. Inspired by [20] we consider the problem when the initial data of (p[, w[, E'^) is 
around the steady state (/>*, 0, 0) and make use of energy estimates to analyze perturbation 
of the global in-time solutions. To this end, we employ the fourth-order wave equations 
for ^/pl and the equation of the vorticity of velocity The poisson equation is used to 
deal with the coupling of the two carriers and some technique is used to deal with the 
smallness both of e and r. 

Since we are interested in not only the global existence theory but also the asymptotical 
analysis of strong solutions with respect to small parameters, we deal with the scaled IVP 
problem fll.lOI) - fll.l2l) and (12.31) directly. Because the scaled scaled IVP problem fll.lOp - 
fll.l2p and (12. 3p is equivalent to the original IVP problem (ll.3p - (ll.5p and (12. ip for strong 
(classical) solutions. For simplicity, we take A = 1 and let {.)t denote dt{.) and omit the 
index e, r to simplify the presentation in the following argument. From (I1.10p - (ll.l2p and 
(12. 3p the equations for i/ji = -y/pf with Ui = ulii = a, h) can be obtained 

r'Au + ^..+^^ + ■ i^^E) - ^V^(^^. u,) 

' AP.m + ^-'-^-O, (3.1) 



A 4V>, 
with the initial value 

?/^i(x,0) := tpi^ix) = ijlix) = \^pjx), 

Mx, 0) := M^) = ■ < - < ■ v^;. 

Also from (I1.10p - (I2.3I) with the fact {ui ■ V)ui = |V(|uip) — x (V x m), the equations 
for Ui = ui (i=a,b) 



T^u^t + u, + -V(|n.|') - T^u, X + = q,E + (3.2) 



where (pi = V x Ui denotes the vorticity of Ui. Taking curl of (3.2), we have 

T^(j)it + + T^iUi ■ V)0, + T^(j)iV ■ U, - T\(t), ■ V)Ui = 0, (3.3) 

Introduce new variables Wi = i/ji — with i = a,b, then the system for {wa, Wb, (pa, 4>b, E) 
is 

e^A^w 1 

T^Watt + Wat + + 2 K + V^)V ■ E - P'M)^Wa = faU (3.4) 

r'^Wbtt + wu + ^^^^ - \{wb + v^)V ■ E - Pi{pl)/\Wb = Ai, (3.5) 

T^^at + 0a = /a2, (3.6) 
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and 
where 



r^0M + (pb = fb2, 



V-E = wl-wl + 2^^Wa-2y^,Wb, V X E = 0, 



fii ■■= fiiix,t) 



-T W, 



q,Vw,E + imw, + vV:)') - ^'(pD) Am 



Wi + VP* 



2(m, + v^o*) 



/i2 := fi2{x, t) =r^((0j ■ V)Mj - {ui ■ - -Ui), i = a, b. 
The last term in (13. 9p can be decomposed by using equation (11.101) as 

=T {-WitV ■ Ui - 2ui ■ Vwit - — ^ + Vwi ■ {{ui ■ V)Ui 



+ ^ 2^ \9kUi\ 

k,l=l 



2(m. + VP*) 
(wi + Vpf) 



— Mi ■ V('Uj ■ V^i 



+ 



1 



Wu + Ui-Vwi){ui-Vwi)}, i = a,b. 



2{w^ + ^ 
The initial conditions for (I3.4p - (I3.7I) are 

Wi{x, 0) := Wio{x) = - a/p*, 0) := 0io(a;) = -V x Mio(x), 

r 

^^4(2;, 0) := mil (a;) = -{-Uio ■ Vwi^, - -(mj,, + ^/p*)V ■ WiJ, i = a,b. 

r Z 

We will also use the relation between V ■ Ui and Vmj, Wu 

2wit + 2Mj ■ Vwi + (mj + a/p^)V ■ = 0, i = a,b. 

3.2 The a-priori estimates 



(3.7) 
(3.8) 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



In this section, we will mainly study the reformulated equations (I3.4p - (I3.8I) in order to 
obtain the a-priori estimates of Wa, Wb, 4>a, 4>bi E. 
Set the workspace as 



X(T) = {{wa, Wb, Ua, Ub) E L°°([0, T]; {H^{R')y x {n'{R')Y} 
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and assume the quantity 

St = max{|| {wa,Wb)i.,t) \\l + || T{dtWa, dtWb)i.,t) \\l + \\ T{ua,Ub){.,t) 11^4} 

+ / {II iUa,Ub)i.,t) 11^3 + II iWa,Wb)i.,t) \\l + || f^^jdt, (3.13) 

Jo 

is small, then by Sobolev embedding theorem we know that the sufficiently small St can 
assure the positivity of ipa^'^^b as 



2 - ' V ra - 2 V ra> ^ - ' V - 2 

By Sobolev embedding theorem, from the assumption for 6t, we also have 

II (D'^wa, D'^Wb, rD^Wat, rD^Wbt) \\L^iR^>,[o,T])< cdT, \a\ < 2, \P\ < 1. (3.14) 
II {TD°'Ua,TD°'ub) \\l°-{r3x[o,t])< c5t, \a\ < 2. (3.15) 

/ II {D^'Ua, D°'Ub,T'^Uat,T'^Ubt){.,t) ||ioo(^3) dt < c6t, \a\ < 1. (3.16) 

Jo 

The last inequality 03.161) is obtained from the equations for Ua, Ub and Sobolev embedding 
theorem, the assumption for 5t- The c or C denote the generic positive constant and does 
not necessarily be the same here and after. Using Lemma 12.71 from the poisson equation 
(13. 8p we have 

II E ||L°°([o,T];7T!5{i?3))< c5t, \\ D"' E \\ (r3 x[o,T])< c5t, \a\ < 3. (3.17) 

Next, we will establish energy estimates to extend the solution to global one. 
We have the main a-priori estimate lemma. 

Lemma 3.2 Suppose {wa,Wb,Ua,Ub, E){x,t) is local solution with 6t <^ 1, then it holds 

Ei{t) + [ E2{s)ds < cAo, (3.18) 
Jo 

for t G (0, T) and c > is a constant independent of e and r. The Aq is defined in 
Theorem 2.6, and here 

Ei{t) :={|| {wa,Wb){.,t) \\1+{t + £^) II {D'wa,D'>Wb){.,t) f 
+Te^ II (DV,^V)(.,t) f II iwat,Wbt)i.,t) \\l 

+ II {D^Wat,D^Wbt)i.,t) f +r' II {Ua,Ub)i.,t) 

+t' II {D'>Ua,D'>Ub){.,t) f + II E{.,t) 11^,,}, 
E2{t) := {II {Vwa,Vwb){.,t) Wl+e^ \\ {D^Wa, D^Wb){.,t) f 

+ II iWauWbt)i.,t) \\l +T II {D^Wat,D^Wbt)i.,t) f 

+ II iUa,Ub)i.,t) 11^4 +r II {D''Ua,D'>Ub)i.,t) f + II E{.,t) II^J. 
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Proof: Step 1. The estimates for Wa,Wb- 
Step 1.1. basic estimates. Assume r < 1 for simplicity. Multiplying (13 ■4p by {wa + '2^Wat) 
and fl3.5p by {wb + 2wbt), integrating by parts the resulted equations over R^, summing 
the resulted two equalities and noticing the facts from equation (13.81) 

= - ! \V-E\'dx-][ V{wl-w!)-Edx, 
4 Jr3 4 Jji3 

and 

/ dii^a + v^) V ■ E)2wat - \{{wb + v^)V • E)2wu}dx 
= ~ [ \V-E\'dx, 

4:dtjR3^ 

then after a tedious but straightforward calculation we have 

,2 ,,,2 



d_ 

It 



f VO ID 

+P;,{pl)\Vwb\^ + ^(|A«;,|2 + \/\wb\^) + i|V ■ E\^}dx 



+ / {(2 - r^){wl + wl) + P:(p:)|V«;.|2 + Pl{pl)\Vwb? + ^(I^^^P + |A«;,| 
Jb? 4 



+^|v-£;p}(ix 

4' 



2 JR3 



+ / {/ai(a;,t)(wa + 2wat) + /6i(2;,t)(^ya + 2i(;fet)}(ia;, (3.19) 

The right-hand side of (13.191) can be analyzed as follows. By Sobolev embedding theorem 
and Holder inequality, Young's inequality 

WiVwi ■ Edx < II Wi \\lA\ I|l2 ■ II E 11^6 

<c(|| w^ 11^2 + II Vwi lUOdl ||l2 ■ II E ll^e) 

<c(5t)^(|| V^/;,f+ II V-i?f), (3.20) 

i=a,b. Here we have used Lemma [2. 71 to estimate || DE p by || V ■ i? ||^. Some other key 
terms of the right-hand side are analyzed as 

/ [Pl{{w. + ^f) - Pm]^w, ■ {2wu)dx < c{5t)--{\\ Aw, f + W Wu f ), (3.21) 
Jb^ 

12 



/ T'ui ■ Wwu{2wu)dx = - r'V ■ Ui^wnYdx < c(5t)^ || w,t \\\ (3.22) 



T UiW{ui-Wwi){2wit)dx < c{5t)^ \\ {Vwi, Awi,Wit, \\ , (3.23) 

In (13:231) we have used || Dui f< c(|| V ■ Wi f + || V x m f), and || Vwi f , || Wu f 
to estimate V ■ Ui through equation fl3.12p . Then, by fl3.19p - fl3.23p . using integration by 
parts, Holder inequahty. Young's inequahty and the Moser type Lemma 12. 9[ Lemma 12.101 
to estimate the other terms of the right-hand side of fl3.19p , we can arrive at 

d f ViP' w'^ 

+ Plipl)\Vw,\' + j{\AwX + \^w,\') + ^|V ■ E\'}dx 
+ [ {(2 - r')iwl, + wl) + P'Apl)\^Wa? + Pi{pl)\^w,\' + + |Au;,n 



ij3 

+ ^|V ■ E\^}dx 



4 



4' 

< c(5t)^ II (VW„ VWi,, Wat, Wat, V ■ E, (j)a, (pb) f +c((5t)^ || (Aw„, AWfe) f , (3.24) 

The right hand side of estimate fl3.24l) will be used later in the closure of the a-priori 
estimates. 

Step 1.2. the higher-order estimates for Wa,Wh. Differentiate 03.40 and fl3.5p with 
respect to x, then the functions Wa := D'^Wa,Wb '■= D°'Wh and E := D"i?(l < |a| < 3)E 
satisfy 

2 

T^Witt + Wit + ^^^Wi + ^{Wi + v^)V ■ E - P',{p*)AWi 



=Dya{x,t) - D^qiwi + v/p*)V ■ E) + |(«;, + ^p*)V ■ E 

clef 

= Fi{x, t),{i = a,b, qa = l,qb = -1.) (3.25) 

Multiplying (13.250 for i = a by {w^ + 2wat), and (13.250 for z = 6 by {wh + 2wlt), integrating 
by parts over R^, summing the resulted equalities, also noticing the facts 

\ {l^{^a + V^)V ■ [E){Wa + 2Wat) " T^(Wf> + aM)^ " (^)(«'6 + 1Wbt)}dx 

^ ^ 'V-E\''dx + \^ [ \\/-E\^dx-\ [ V ■ED''{wl~wl)dx 



4 .7^3 4 dt JR3 ^ Jji3 

If V-ED''{wl-wl)tdx+j ]-WaV ■{E){^a + 2Wat)dx, 



* we can first assume the solution (wa,Wb,Ua,ui,) has higher order regularity so that we can take 
derivatives since the final a-priori estimation will be still valid for these solutions by applying the Fricdrich 
mollifier to {wa,Wb,Ua,Ub) ■ 
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IwbV ■ {E){Wb + 2Wbt)dx, (3.26) 
after a tedious but straightforward computation one can get 

^ / {'^'^w'at + r'^uTauTat + \w'a + r'^Wb] + r'^WbWbt + ^wJJ^ + (p* ) | VuJ; | ^ 

+ n'(p^)|V^Z^p + j(|AS};|2 + \/\Wb\^) + i|V ■ E\^}dx 

+ / {(2-r2)(^^ + 5y^^)+P^(p:)|V.?;|2 + A'(pD|V5^r 

1 ~ 

+ — (|A{?;|2 + |A?^n + -|V ■ 



{Fa ■ {wa + 255;0 + Fb-{Wb + 2Wbt)}dx + - / V • - wi)dx 

2 Jr? 2 
+ ^ / w;feV-^(S}; + 2S5;i)rfa;, (3.27) 

2 Jij3 



Similar with the analysis of basic estimates, using Moser type inequality Lemma 12.91 
Lemma [2. 101 and the priori assumptions fl3.13p - fl3.17p and using Holder inequality, Young's 
inequality to estimate the terms of the right-hand side of fl3.27p . we can arrive at 

^ j {r^Wa] + r^WaWat + ^Wa^ + T^Wb] + ^'^WbWbt + + ^a(Pa)|V5^P 

+ Pi{pl)\VWb\' + j(|A5};|2 + \/\wb\^) + i|V ■ E\^}dx 

+ / {{2-T^){Wa] + Wbl) + P'a{pl)\VWa\^ + P;,{pl)\VWb\^ 

e"^ 1 ~ 

+ j(|A{?;|2 + \A^b\^) + -|V • E\^}dx 

< c4 II [VWa, VWb, Wat, Wbt, <Pa, (pb, ^ ' E) \\l +c4 || {D^Wa, D^Wb) f, (3.28) 



Note that, we can not deal with the last term in (13.281) by the energy of left-hand side 
now, so we have to do the highest-order estimates in different way in order to overcome 
the difficulty. 

Step 1.3. the highest-order estimates for Wa,Wb. Taking |a| = 4, we can get the 
equations for u?^ := D°'Wa,Wb '■= D'^Wb and E := D"E. We also use the form of (13.251) 
for simplicity. This time, using {wa + 2TWat) to multiply fl3.25p ,-^ and (wl + 2TWbt) to 
multiply fl3.25l) j->, but for |a| =4. we can get as former 



d 
It 



I {r^Wat + r'^WaWat + ^Wa^ + r^Wbf + r'^WbWbt + ]:Wb^ + TP^{pl)\VWa\ 

Jr3 2 z 
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+ rPl{pl)\Vw,\' + ^{\Awa\' + \Aw,\') + ^1 V ■ E\'}dx 

+ 'l(\AwX + \Am,\') + ^\V-E\'}dx 
[ {Fa ■ i^ITa + ^rWat) + Ffe ■ (5}^ + 2TWu))dx I V ■ ED^'iwl - wl)dx 
+ 1 / ■ ED^'iwl - wDtdx -I WaV ■ E{Wa + 2TWat)dx 

^ JR^ ^ Jr3 



'i?3 ^ JR3 

+ 1 I WbV ■ E{Wb + 2TWbt)dx, (3.29) 

2 .Ir3 



In the right-hand side of (13.291) . the terms multiphed by 2TWat,2TWbt need a special 
analysis. Taking i = a for example, the key terms are analyzed as 



'i?3 

^ r 1 

L-TaUU^a ^ y Pa) ) - J^a\Pa)\ \ ^ '^a\ "T "^"t II ^ "^a l' 

+ c4r II ^D;if, (3.30) 



/ r'^UaVWat ■ 2rWatdx = - T^V ■ Ua\Wat\'^dx < cd^T \\ Wat (3.31) 

Jr3 Jr3 



and 



r n„V(^Za ■ Vii;a) ■ 2TWatdx 

RJi 



<--T: I rijUa-WWafdx + I 2T^{Ua-VWa)UatVWadx 

at jfjji jjis 

+ c4 II VWa f +c4^ II Wat f , (3.32) 

The other terms in the right-hand side of f l3.29p can be analyzed just use Moser type 
Lemma [2^ Lemma [2. 101 the assumptions f l3.13p -( l3rT7l) and the Sobolev embedding theo- 
rem, the Holder inequality. Young's inequality. In a word, these estimates with the above 
estimates (l3:29|) - (13321) will lead to 

^ f f 3 2 I 2^ I 2 I 3- 2 I 2' ■ , ^ 2 , 7-)// *Mv7' 12 

— / {t Wat+r WaWat + -Wa + T Wbt + T WbWbt + -Wb + T P^ipJlV Wa\ 

at J Z Z 



TPl{pl)\VWb\' + ^{\Aw,\' + \Awb\') + ^|V • E\'}dx 



+ I r[P:,{{Wa + ^af) - P'aipim'^a^dx + jj^^ r{rUa ■ VWafdx 
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< c4 II (Vw;^, ViOfo) 11^ +c4£^ II f +c4 II V ■ ^ 11^ 

+ c4r II (D^Wat, D^Wbt) ir +c4 II (0a, 0b) II4, (3.33) 

Note that the right hand side of the estimates (13.241) . fl3.28l) . fl3.33p will be treated later in 
terms of the estimates of (j)a,4>b- 

Step 2. The estimates for (pa, <pb- Differentiating the equations (13. 6p and (13. 7p for 
(pa,<Pb with respect to x, then (pa = D°'(pa,(pb = -D"06(|«| < 4) will satisfy (taking (pa for 
example) 

T^^at + 0a = I^°/a2, (3.34) 

recall fa2 in (I3.10p for i = a. Taking inner product between 2(pa and (13.340 . integrating 
over R^, we obtain 

r'f^f \^a\'dx + 2 [ \^a\'dx= [ D''fa2-2^adx, (3.35) 



R3 JR3 JR3 



The terms in right-hand side of (I3.35P can be estimated using Moser type Lemma | 
Lemma 12.101 Young's inequality and the assumptions (I3.13p - (l3.17p and the inequality 
II Du \\< c(|| V ■ M II + II V X w II) and also the presentation of V ■ by Wat, Vtfa through 
equation 13.121 for i = a. Then we deduce 

d f ~ f ~ - - - 

'^'^Zu / I'Pal^dx + 2 \(pa\'^dx < C6^ II (pa \\l +c6^T \\ Wat \\l +C(5| || Vl^a II4 (3.36) 

at Jj^s J ^3 

Step 3. The closure of energy estimates. The assumption 5^ -C 1 and the 

combination of the estimates (13.240 . (13.280 for all |a| < 3, and (13.330 for all |a| =4, (13.360 
for all |a| < 4 can give us 

jHi{t) + H2{t) <Y,t\\ Uii.,t) |Uoo(^3) ■ II T\ui;t) |Uoo(K3) ■ II D'wi f (3.37) 

i=a,b 

where Hi{t), H2{t) are two terms satisfying 

< ciEi{t) < Hi{t) < C2^i(t); < C3E2it) < H2{t) < 04^2 (t). 

for t G [0,r], and Ci, 02,03,04 are positive constants independent of e, r, the Ei{t),E2(t) 
are the terms defined in the Lemma [3.21 From (I3.37P we can write 

j^H,{t) + H2{t) < cg{t)H,{t), t e [0,r]. (3.38) 
16 



with 

9{t) = II Ui{-,t) |Uoo(^3) ■ II T'^Uit{.,t) \\L°-iR3) 

i=a,b 

The assumption fl3.13p then fl3.16p with the Gronwall inequahty applying to 03.381] makes 
us know 

Hi{t) < ce^o9Wrf^ff^(o) < ce^^^HiiO) < CHi{0). (3.39) 
for t G [0,T] provided 6t <^ 1- Integrating fl3.38p on [0,t] and using (13.390 . we derive 

[ H2{s)ds < Hi{0) + Hi{t) + C6tHi{0) < C'Hi{0), (3.40) 
Jo 

The above constants c, C and C denote the positive constant independent of the param- 
eters £ > 0, r > 0. 

It follows from (I3.39l) . (l3.40p and the equivalence between Hi{t) and Ei{t), and between 
H2{t) and E2{t) the conclusion stated in Lemma |3.2[ Thus the proof of Lemma 13.21 is 
completed. 



3.3 The global existence and asymptotical limits 

The proof of global existence (Theorem \2.1\) : The Theorem 12. II is a direct conclusion 
of the combination of the local existence theory Lemma 13.11 and global a-priori estimates 
Lemma 13.21 in terms of the variable transformation presented above and the standard 
continuity argument, we omit the details. 

The proof of semiclassical limit (Theorem \2.3\) : Starting from Lemma 13. 2[ using a 
continuity argument, one can easily prove the existence of the global in-time solutions of 
the original problem ( ]1.10p -( !273l) with any small e and r provided the Ai > then Aq > 
small enough. 

Let {ipl,'i/jl,ul,ul, E^) be the solution of fll.l0p - fl2.3l) . then from Lemma \^72\ and the 
poisson equation (13.80 the uniform estimates to e hold 

1 1 

E II i^tira- ^M)^d^irb- lit. +E n 

k=0 k=0 

+ II E%.,t) 11^6 < cAo, (3.41) 



{II m - VPI)^ i^l - VpI)){;S) U + W {dtra,dtrt){;S) \\i}ds < cK,t (3.42) 

/ {E II (5M,5M)(-,^) 11^.- +E II i^tE^i-.s) \\l.-.}ds < cAo (3.43) 

•^^ k=0 k=0 



for any t > 0. The right-hand sides of the above inequalities are independent of e. Thus 
these uniform estimates and Aubin's lemma imply the existence of subsequence denoted 
also by {i/jljipl^uljUl, E"^) such that 

i^l^^a, ^I'l^A m C{0,t;C!nHf-^{R')), (3.44) 
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ul^Ua, ul^Ub in C{Q,t-Clnn'l;:{R')), (3.45) 
E^^E in C{0,t;Ctnnlr{R')), (3.46) 

with s G (0, |), as £ ^ 0. We also have 

yV(^)-^0 m L\0,t;HUR')) 

as e 0. Thus fl3.4ip - fl3.46l) allow the e pass to the zero, and the limiting solutions 
satisfy 

2tPadt^a + V ■ i^PlUa) = 0, 

r'dti^lUa) + T^Vi^lUa ® Ua) + Pai^l) + - ^IE = 0, 

2^,dtiJb + V ■ {iPlu,) = 0, 
T^^t{^lu^,) + T^V{^lu^, ® u^,) + VPfe(^b') + ^/^fe^Wfe + ^lE = 0, 

Let Pa = ("^a)^, Pb = K is easily to verify that {pa, Pb, E) solves the bipolar hydrody- 

namic model fll.6p - fll.8p . The convergence of the bipolar quantum QHD model to bipolar 
hydrodynamic model is established, and the proof of the Theorem 12.31 is complete. 

The proof of combined semiclassical and relaxation limits (Theorem \2.4\) : Since 
the estimates established for the solutions in Lemma 13.21 hold uniformly for any small e 
and r, thus we can study the combined limits as both e and r tends to zero freely. 
Let (4"''\#'^ 

,u'"J''^\ul^'^\ E'^'^''^'^) be the global solution derived in Theorem 2.6, by 
the estimates f l3.18p in Lemma 13.21 we have the uniform estimates about e and r as 

II (^i^'^^ - v^,#'^ - v^)(.,t) 11^+11 (™^^),™^^^)(.,t) ||^4< cAo (3.47) 
II (ra,^^^),ra,V^^^))(.,t) 11^+11 E(^'^)(.,t) ||?,.< cAo, (3.48) 

and 



II (^^^^ - ^/Pl^i^t^ - V^)(-,^) lis + II i^,^p!^''\^,i;i^''^)i.,s) \\l)ds 
<cAot, (3.49) 

All iul^^'\u^;^'^)i.,s) 11^4 + II E^^'^\.,s) \\l,.)ds < cAo, (3.50) 
Jo 

for any t > 0. 

Also use Aubin's lemma with the above uniform estimates, we can get the subse- 
quence(not relabelled) and functions denoted also by ■0a, ipb, Ua, Ub, E such that as e, r — > 

V^^^) ^ ija, 4^'^'^ - i^b ^n C{0,t; Cl H Hf-^R^)), (3.51) 
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u^a'^'^^Ua, v!-^'''^ Weakly in L'^{0,t;n^{R^)), 

Ei-'^) m C(0, t- Cl n 



(3.52) 
(3.53) 



for any t > and s G (0, \). 




0. 



(3.54) 



Thus the above converging results allow the solutions to pass to the limit r, £ — from 
the bipolar QHD model to the bipolar drift-diffusion(DD) model: 



which is equivalent to the bipolar DD model fll.15p - fll.16p in Section 1 for strong solution. 
Namely, (pa = (V'a)^, Pb = (V'b)^; E) solves the bipolar Drift-Diffusion model f ]1.15p -f fLT6l) . 
The proof of Theorem 2.7 is completed now. 
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